Circularly polarized modes in magnetized spin plasmas 
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Abstract 

The influence of the intrinsic spin of electrons on the propagation of circularly polarized waves in 
a magnetized plasma is considered. New eigenmodes are identified, one of which propagates below 
the electron cyclotron frequency, one above the spin-precession frequency, and another close to the 
spin-precession frequency. The latter corresponds to the spin modes in ferromagnets under certain 
conditions. In the nonrelativistic motion of electrons, the spin effects become noticeable even when 
the external magnetic field Bq is below the quantum critical magnetic field strength, i.e., Bq < 
Bq = 4.4138 X 10^ T and the electron density satisfies uq ^ Uc — lO^^m"'^. The importance of 
electron spin (paramagnetic) resonance (ESR) for plasma diagnostics is discussed. 
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I. INTRODUCTION 



During recent years there has been a rapid increase in the interest of quantum plasmas, see 
e.g., Refs. [THIS]. This has been stimulated by experimental progress in nano-scale plasmas 
[0], ultracold plasmas [IB], spintronics [T7j, and plasmonics [18j. However, already more than 
forty years ago, lannuzzi pTS] established the possibility for the existence of electron spin 
(paramagnetic) resonance (ESR) in a fully ionized low-temperature plasma, and predicted 
its importance, e.g., in the plasma diagnostics for measuring the particle density with a 
greater precession than the conventional technique, in determining the particle velocity 
spectrum perpendicular to the magnetic field, as well as to calculate the magnetic field in 
the propagation of electromagnetic (EM) waves (e.g., whistlers, Alfven waves, shock waves 
etc.) in plasmas. Recent investigation [SD] along this line indicates that besides the currently 
prevalent laser methods, ESR technique can successfully be used for plasma diagnostics, e.g., 
measuring the electron densities in the microwave region. Furthermore, the importance of 
spin effects in plasmas has also been studied using kinetic plasma theory [TlHTSt [2n 122] , 
with applications to wave propagation p!THT3] as well as other phenomena [211 [22] . The 
hydrodynamic description of spin plasma waves can also be found in the literature (see, e.g. 
Refs. [M El 123]). 

Although, whistler waves have been studied for almost a century, they are still a subject 
of intense research in view of its importance not only in space plasmas, but also in many 
astrophysical environments, e.g., in the atmosphere of neutron star envelope, in the coherent 
radio emission in pulsar magnetosphere etc. Quasilinear theory [21] and simulation [25] show 
that whistler waves can be used to resonantly accelerate electrons. Furthermore, they can 
also be used to interpret the fine structure of Zebra-type patterns and fiber bursts in solar 
type II and IV radio bursts [26j. Thus, the occurrence of ESR might be useful for the electron 
acceleration in the propagation of EM radiation in the pulsar magnetosphere as well as for 
plasma diagnostics in the microwave region in laboratory experiments if the conditions favor. 

In the present work, we will derive and analyze the dispersion relation for the propagation 
of circularly polarized (CP) EM (CPEM) waves in a magnetized spin plasma using a spin 
fluid model. Various fluid models are appropriate in different regimes (see e.g., Refs. [QllTO]). 
The model to be used contains the Bohm-de Broglie potential, the magnetic dipole force and 
includes the spin-precession dynamics as well as the spin magnetization current. Its basis 
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can be found in e.g., Ref. [10], and more rigorous foundation can be given starting from the 
kinetic theory presented in Ref. [13]. Specifically, we will focus our discussion to the ESR 
as well as the properties of spin modified whistler-like modes. 

II. WEAKLY NONLINEAR EVOLUTION 

The nonrelativistic evolution of spin— 1/2 electrons can be described by [10] 

e fV'^ /rr\ / 2u \ 

+ v.. V) V. ^ -- (E + V. X B) - VPjm.n. + —V [-^) + j S. VB, 

(1) 

(ai + v,.V)S = -(2;u/n)(BxS), (2) 

where ne, me, Ve respectively represent the number density, mass and velocity of electrons, 
E (B) is the electric (magnetic) field, Pg is the electron pressure. Also, S is the spin angular 
momentum with its absolute value |S| = h/2\ fi = — {g/2) /i^, where g ~ 2.0023193 is the 
electron (^-factor and fiB = eh/2me is the Bohr magneton. The equations are then closed 
by the Maxwell equations. 

V X E = -dtB, V.B = 0, (3) 

V X B = /io {6odtE-eneVe+ V x n^S) , (4) 

The above equations ([T]) and ^ apply when different spin states (i.e., spin- up and spin-down 
relative to the magnetic field) can be well represented by a macroscopic average. This may 
occur for very strong magnetic fields (or a very low temperature), where generally the lowest 
energy spin state is populated. Alternatively, when the dynamics on a time-scale longer 
than the spin-flip frequency is considered, the macroscopic spin state is well-described by 
the thermodynamic equilibrium spin configuration, and the above model can still be applied. 
In the later case, the macroscopic spin state will be attenuated by a (thermodynamic) factor 
decreasing the effective value of |S| below h/2. However, this case will be not considered 
further in the present paper. As a consequence, our studies will be focused on the regime of 
strong magnetic fields, as seen in astrophysical plasmas [27] . 

In what follows, we will assume the propagation of a CPEM waves to be of the form E = 
(x ± iy) E{z, t) exp{ikz — iujt)+ c.c, along an external magnetic field B = Bqz, where E{z, t) 
is the weakly modulated wave amplitude (i.e. we assume | (1//) df/dz\ <^ k,\{l/ f) df/dt\ <^ 
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for all variables /), k (u) represents the EM wave number (frequency) and c.c. denotes 
the complex conjugate. In the interaction of high-frequency (hf) EM waves with the hf 
electron plasma response, the use of cold plasma approximation is well justified in view 
of the fact that for large field intensities and moderate electron temperature, the directed 
velocity of electrons in the hf fields is much larger than the random thermal speed. Moreover, 
it can also be shown that the density perturbation associated with the high-frequency (hf) 
EM wave is zero. Thus, taking the curl of Eq. Q and using Eqs. ([2])-(|4]) we readily obtain 
the following evolution equation for CPEM waves. 



= —dtB + —dt ( a^B + — V^e X a^E I - v^,V x 9zve + —dt 
me eue \ rie J e/io 
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+ — ^ — V X [(V X B) X B] + X {S'^VBa) 
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^° -V X {dt^ X B) X [(V X rieS) x B] (5) 



meUe mehUe 

The weakly nonlinear equation ([s]) is a useful result when considering the interaction between 
low-frequency (If) and hf fields, where the If fields are induced by the ponderomotive force. 
Equation ([s]) then needs to be complemented by equations for the If dynamics, and naturally 
the number of dependent variables can be further reduced. However, before this line of 
research is pursued, a more through study of the linear theory should be made, as will be 
undertaken in the next section. 



III. LINEAR THEORY OF WHISTLER WAVES 

Introducing the variables B± = i: iBy, E± = i: iEy etc., suitable for CP waves, 
and limiting ourselves to the linearized theory, we obtain respectively from the Faraday's 
law and the spin-evolution equation [15]. 

B± = ±—E±, S± = Tt—— — r (6) 
00 n[uj ± ujg) 

Using Eq. ^ to express the free current as well as the magnetization current in terms of 

the electric field, and using Eq. Q, the following linear dispersion relation is obtained for 

the CPEM modes 

4=1- < - f^yi^'i 
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which can be rewritten as 

nl(l + ^^)=l- (8) 

where = ck/u is the refractive index, and where the upper and lower sign respectively 
stand for the left-hand circularly polarized (LCP) and right-hand circularly polarized (RCP) 
waves. Also, = g'^h/SrrieXl is a frequency which involves the spin correction due to plasma 
magnetization current and Ag = c/upe is the electron skin-depth (inertial length scale). 
Moreover, Upe = ^/noe'^ /eorrie, coc = eBo/rrie and Ug = {g/2)uJc are respectively the electron 
plasma, cyclotron and the spin-precession frequency. In absence of the spin-magnetization, 
the well-known classical dispersion relation, namely u"^ = c^k"^ + uJUp^/ {u ± Uc) is recovered. 
A few comments are in order. The first and the second term in the right-hand side of Eq. ([T]) 
appear due to the displacement current and the free electron current. The term involving 

appears even in absence of the external magnetic field, since the spin perturbation is due 
to the wave magnetic field not the constant field Bq which does not provide any magnetic 
dipole force. Note, however, that an unperturbed spin state with |So| = h/2, that has been 
used in the derivation, typically requires that the external magnetic field is strong. Thus, 
inclusion of the electron spin perturbation leads to a modification of the dispersion relation 
for transverse plasma oscillations. This modification is clearly substantial when u < Ug ^ 
w^, i.e., when hiOc > mgC^ for tUpe < V2ujc, where c is the speed of light in vacuum. This 
corresponds to a regime of very strong magnetic field in which the external field strength 
approaches or exceeds the quantum critical magnetic field, i.e., Bq > Bq = 4.4138 x lO^T. 
In such a situation relativistic effects might be important. On the other hand, for the 
nonrelativistic motion of electrons we have hwc < meC^, i.e., Bq < Bq for cOpe > a/2 cUc- In 
this case, the density regime in which the magnetic field is 'nonquantizing' and does not 
affect the thermodynamic properties of the electron gas, is no ^ Uc — lO^^m"^ and the 
temperature Tg ^ Tb — huc/kB, where fcs is the Boltzmann constant. Thus, in the strong 
magnetic field and highly dense plasmas, the electron spin effect can no longer be neglected, 
rather it modifies the wave dispersion leading to new eigenmodes. 

Inspecting now the term oc h in Eq. ([T]), we note that 

So, the spin current can be much larger than the classical free current when \Jm±I Jc±\ ~ 
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hk"^ /iTLeUJ ^ l.This basically holds when either (i) hujc ^ fneC^, u < Uc, ck or (ii) huc < meC?, 
uj < ujc and u <^ ck is satisfied. The case oi u > Ug > Uc is rather less important, as it 
does not give rise to wave propagation (ri^ < 0). Also, in the very If regime u <^ Wc, the 
ion motion can be of importance, and we will therefore not consider that case. Thus, one 
important mode could be the RCP If {uj < Uc) EM waves (whistlers). Let us now see how 
the dispersion relation reduces when either of the two cases is considered. In the limit of 
\Jm±/ Jc± \ ~ hk'^/meOO ^ 1, the dispersion equation (|8| reduces to 



from which one finds for u ck a purely spin-modified frequency Ui ~ l^^^g — Also, 
if <^ ck and ^ Ug, we have U2 ~ T^g and ~ c^A;^. The frequencies 001^2 may 
correspond to the spin waves in ferromagnets under certain conditions j29j. Figure 1 dis- 
plays the modes for RCP waves obtained as numerical solutions of the dispersion equation. 
Evidently, there exist two eigenmodes apart from a hf (w > Ug) one, namely a mode close to 
the electron-cyclotron or spin-precession frequency, and the other one is the If mode below 
the cyclotron frequency. In contrast to the hf modes, the spin modified If modes propagate 
with the frequency below that in the classical case. On the other hand, in the very hf regime 
{uj ^ Ug), ^ 1 > 0, and so we have u"^ = c'^k'^, which is the dispersion relation of a 
EM wave in vacuum. Evidently, since n'j^ < for uj > Ug, there must exist an intermediate 
frequency uj{> ujg) at which the solution for must pass through a zero value, and becomes 
positive again. Thus, the cut-off frequencies for which n|. = are obtained as 



where =F stand for RCP and LCP waves respectively. Clearly, the RCP waves have lower 
cut-off frequency than the LCP modes. On the other hand, the resonances for the RCP 
waves ( n\ — )■ 00) associated with both the orbital and the spin-gyration, occur (LCP mode 
has no resonance) as either uj ^ uJc 01 uj ^ Ug, i.e., when the angular frequency of the wave 
electric field matches either due to electron cyclotron motion (cyclotron resonance) or due to 
the intrinsic spin of electrons (ESR). At the resonance, the transverse field associated with 
the RCP wave rotates at the same velocity as electrons gyrate around Bq. The electrons thus 
experiences a continuous acceleration from the wave electric field, which tends to increase 



{U ± Ug) (1 



u^/c^P) +ujf, = 



(10) 
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FIG. 1. Different eigen modes for RCP waves obtained as numerical solutions of the dispersion 
equation ([T]) are shown with respect to the normalized wave number and frequency for Bq = 
5 X 10^ T < Sq, no = 7 X 10^^ m-^ > n^. 

their perpendicular energy. Therefore, it is not surprising that RCP waves propagating along 
the the external magnetic field and oscillating at the cyclotron frequency or spin-precession 
frequency are absorbed by electrons. This may be the consequence to the recently developed 
experiment based on microwave absorption and the ESR to be successfully used for plasma 
diagnostics [20]. On the other hand, since the spin effect is appreciable in the strongly 
magnetized dense plasmas, the ESR could well be relevant for the coherent EM radiation 
in pulsar magnetosphere or magnetized white dwarfs. Let us now see how the group speed 
[vg = duj/dk) of the CP waves is modified with the spin correction. We obtain from the 
dispersion relation [Eq. Q] 



Vg = A/ {2uj/ujI + r) , (11) 



where A and P are given by 



A = -^ + ,^%, r= ^ / ^ \ °l (12) 



2g 2me (W ± UJg) ' (uj ± UJc) 4me {iO ± UJgj 

Clearly, the If {u < Uc) component of a pulse (whistlers) propagates more slowly than the hf 
{u > Ug) components as is evident from Fig. 2. It follows that by the time a pulse returns to 
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FIG. 2. The normalized group speed given by Eq. (11) is plotted as a function of the normalized 



wave frequency for the RCP waves with the same set of parameters as in Fig. 1. 

a ground level it has been stretched out temporarily, because the hf component of the pulse 
arrives slightly before the If components. This also accounts for the characteristic whistling 
down-effect observed at ground level. Moreover, the group speed Vg of the whistlers increases 
in the frequency regime < u < Uc/2, and decreases in the other subinterval Uc/2 < u < Uc 
before it reaches the maximum nearer the resonance point. However, the group speed of the 
hf modes approaches the speed of light as u (> Ug) increases and gets saturated at large u. 
From Fig. 2 we also note that the spin force reduces the group speed in strongly magnetized 
dense plasmas. 



IV. SUMMARY AND DISCUSSION 

To summarize, the dispersion relation for the propagation of hf CPEM waves is obtained 
in a magnetized spin plasma. The electron spin modifies the plasma current density and 
thereby introduces a correction term in the dispersion relation, which, in turn, gives rise a 
new CP hf mode. The spin effects are seen to be substantial in the very strong magnetic 
field {Bo > Bq = 4.4138 x lO^T) and highly dense plasmas (no ^ ric — lO^'^m^^) where the 
relativistic effects might be important. However, in nonrelativistic plasmas, the spin of elec- 
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trons can also be important in the case Bq < Bq together with hq ^ ric. In particular, when 
the spin current dominates over the classical free current the RCP EM waves resonantly 
interact with the electrons only at the spin-precession frequency. Such resonance should be 
helpful for particle acceleration in the coherent radio emission of the pulsar magnetosphere 
or magnetized white dwarfs. The study of the spin modified modes might also be important 
at least from the diagnostic points of view, since the observation of the propagation charac- 
teristics of the wave modes may be used in order to determine the physical parameters in 
plasmas [3D]. Lastly, the electron spin-resonance (ESR) could be an important consequence 
to the recently developed experiment for plasma diagnostics in the microwave region, if the 
conditions favor [20j . 
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